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We compare two classes of hybrid equations of state with a hadron-to-quark mat- 

m : 

ter phase transition in their application to core collapse supernova simulations. The 
first one uses the quark bag model and describes the transition to three-flavor quark 
matter at low critical densities. The second one employs a Polyakov-loop extended 

CO 

Nambu-Jona-Lasinio (PNJL) model with parameters describing a phase transition 
to two- flavor quark matter at higher critical densities. These models possess a dis- 

> : 

J^j , tinctly different temperature dependence of their transition densities which turns 

out to be crucial for the possible appearance of quark matter in supernova cores. 
During the early post bounce accretion phase quark matter is found only if the 
phase transition takes place at sufficiently low densities as in the study based on the 
bag model. The increase critical density with increasing temperature, as obtained 
for our PNJL parametrization, prevents the formation of quark matter. The fur- 
ther evolution of the core collapse supernova as obtained applying the quark bag 
model leads to a structural reconfiguration of the central protoneutron star where, 
in addition to a massive pure quark matter core, a strong hydrodynamic shock wave 
forms and a second neutrino burst is released during the shock propagation across 



2 



the neutrinospheres. We discuss the severe constraints in the freedom of choice of 
quark matter models and their parametrization due to the recently observed 2 M 
pulsar and their implications for further studies of core collapse supernovae in the 
QCD phase diagram. 



I. INTRODUCTION 

Stars more massive than 8 M explode as core collapse supernovae, with kinetic explosion 
energies of the ejected material on the order of 10 51 erg. The remnants, the protoneutron 
stars (PNSs), are initially hot and lepton-rich and cool via deleptonization during the first 
30 seconds after the onset of the explosion. Above a certain progenitor mass threshold on 
the order of 40 M Q , which is an active subject of research, stars will no longer explode. 
Such models will collapse and from black holes. The critical mass for a PNS to collapse and 
from a black hole is given by the equation of state (EoS). The commonly used EoS in core 
collapse supernova simulations are based on pure hadronic descriptions, e.g. the compressible 
liquid-drop model with surface effects l| and relativistic mean field theory including the 
Thomas- Fermi approximation for heavy nuclei [2} . The conditions that are obtained in PNS 
interiors during the simulation, are several times nuclear matter density, temperatures on 
the order of several tens of MeV and a low proton-to-baryon ratio given by the electron 
fraction of 1 Y e = Y p < 0.3. Under such conditions, the quark-hadron phase transition is not 
unlikely to take place and the assumption of pure hadronic matter becomes questionable. 
Ab initio calculations of the phase diagram and EoS of quantum chromodynamics (QCD) as 
the fundamental theory of strongly interacting matter come from simulations of this gauge 
theory on the lattice but are still restricted to low baryon densities (chemical potentials). 
In this region of the phase diagram they predict a crossover transition with a pseudocritical 
temperature for chiral symmetry restoration and deconfinement at T c ~ 150 — 170 MeV 
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3, 4]. A critical endpoint for first order transitions is conjectured but lies, if it exists at all, 
outside the region presently accessible by lattice QCD. An interesting conjecture supported 
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1 The proton-to-baryon ratio, Y p = n p /n,B, is equal to the electron fraction, Y e := Y e ~ — Y e + , in the absence 
of muons. 
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by a statistical model analysis of hadron production in heavy-ion collisions and by the large- 
N c limit of QCD suggests the existence of a triple point in the phase diagram [5| due to 
a third, "quarkyonic" phase at temperatures T < T c and high baryon densities p]. This 
state of matter might become accessible in experiments at the planned third generation of 
leavy-ion collision facilities FAIR in Darmstadt (Germany) and NICA in Dubna (Russia) 
3] which thus allow systematic laboratory studies of conditions like in supernova collapse 
and protoneutron star evolution |8[. 

In order to investigate the appearance of quark matter in core collapse supernova mod- 
els, the implementation of a quark-hadron hybrid EoS is required. It must be valid for a 
large range of densities jib, temperatures T and proton-to-baryon ratios. At large baryon 
densities, where lattice QCD cannot be applied due to current conceptional limitations, 
phenomenological models are commonly used. In the present study we will discuss hybrid 
EoS which employ quark matter models that are representative examples from two wide 
classes: bag models and NJL-type models. The popular and simple thermodynamical bag 
model is inspired by the success of the vacuum MIT bag model for the hadron spectrum 
9j. It describes quarks as non-interacting fermions of a constant mass confined by an ex- 
ternal "bag" pressure B. NJL-type models are constructed to obey basic symmetries of the 
QCD Lagrangian like the chiral symmetry of light quarks, and to describe its dynamical 
breaking which results in medium-dependent masses (see [l0(] and references therein). The 
inclusion of diquark interaction channels leads to a rich phase structure at low temperatures 
and high densities with color superconductivity (diquark condensation) in two-flavor (2SC) 
and three-flavor (CFL) quark matter (see, e.g., 11 ] for the phase structure and [12] for the 
relevance to protoneutron stars: detailed information concerning color superconductivity is 



found in review articles 



10, 



13|). These models have no confining interaction and would 



therefore lead to the unphysical dominance of thermal quark excitations already at temper- 
atures well below T c . Their coupling to the Polyakov loop potential is essential to suppress 
the unphysical degrees of freedom and it can be adjusted to fit the behavior of lattice QCD 
thermodynamics at low densities 14j • Extending this effective model to high densities leads 
to the class of PNJL models of which we apply one for this study. 

Different assumptions made for the description of quark matter lead to different critical 
conditions for the onset of deconfinement, which are given in terms of a critical density 
that depends on the temperature and the proton-to-baryon ratio. The resulting different 
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phase diagrams may lead to a different evolutionary behavior for core collapse supernovae. 
The central supernova conditions start from low densities (6 x fm _3 /10 10 g cm -3 ) and 
temperatures (0.6 MeV), where pure hadronic matter dominates, and approach densities 
above nuclear saturation density (0.16 fm _1 /2.7 x 10 14 g cm -3 ) at temperatures on the or- 
der to tens of MeV and Y p ~ 0.2 — 0.3. The nature of the QCD transition is not precisely 
understood on a microscopic level. Therefore, one usually splices independent nuclear and 
quark matter EoS by constructing the phase transition using more or less appropriate con- 
ditions for the phase equilibrium. Popular examples are the Maxwell, Gibbs or Glendenning 
[15I ] constructions, whereby usually the Maxwell construction leads to the smallest region 
between critical densities for the onset and the end of the mixed phase. We will disregard 
finite size effects (pasta structures) and also non-equilibrium effects due to the nucleation 
of the new phase with the justification that weak processes establishing the chemical equi- 
ibrium are fast compared to the typical timescales encountered in supernova simulations 
l^ . We will compare the evolution of a representative core collapse supernova simulation 
in the two different phase diagrams based on the bag model and the PNJL model. The core 
collapse supernova model is based on general relativistic radiation hydrodynamics and three 



flavor Boltzmann neutrino transport 



17H22| 



The manuscript is organized as follows: In Sect. 2 we introduce the standard core collapse 
supernova phenomenology. In Sect. 3 we briefly introduce the two quark-hadron hybrid EoS, 
the quark bag model and the PNJL model. The evolution of a representative 15 M core 
collapse supernova model in the phase diagram is illustrated in Sect. 4 by comparing the 
quark bag and PNJL models. We close with a summary in Sect. 5. 



II. THE STANDARD SCENARIO OF CORE COLLAPSE SUPERNOVAE 

Si-burning produces Fe-cores at the final phase of stellar evolution of massive stars. These 
Fe-cores start to contract due to the photodisintegration of heavy elements and electron 
captures. The latter reduces the dominant pressure of the degenerate electron gas. During 
the collapse, density and temperature rise and hence electron captures, which deleptonize 
the central core, increase. The collapse accelerates until neutrino trapping densities, on the 
order of p ~ 10 n -10 13 g/cm 3 , are obtained after which the collapse proceeds adiabatically. 
At nuclear matter density, the repulsive nuclear interaction stiffens the EoS significantly and 
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the collapse halts. The core bounces back, where a sound wave forms, which steepens into a 
shock wave. The shock wave propagates outwards where the dissociation of infalling heavy 
nuclei causes an energy loss of about 8 MeV per baryon. Furthermore, during the shock 
propagation across the neutrinospheres, which are the neutrino energy and flavor dependent 
spheres of last scattering, additional electron captures release a burst of v e that carries away 
4-5 x 10 53 erg/s on a short timescale of 5-20 ms post bounce (depending on the progenitor 
model and the EoS). Both sources of energy loss turn the dynamic shock into a standing 
accretion shock, already at about 5 ms post bounce. 

The post bounce evolution is determined by mass accretion, due to the continuously 
infalling material from the outer layers of the Fe-core as well as the surrounding Si-layer 
(depending on the progenitor model). On a timescale on the order of 100 ms up to seconds, 
the central density and temperature increase continuously. In order to achieve an explo- 
sion, energy needs to be deposit behind the standing accretion shock which subsequently 
revives the standing accretion shock. Several mechanisms have been suggested, including the 
magnetically- driven 23l-l25|. the acoustic [26] and the neutrino-driven 271 ] . The standard 
scenario, delayed explosions due to neutrino heating, have been shown to work in spher- 



ical symmetry 



28 



29] for the low mass 8.8 M ONeMg-core 30M32l] . For more massive 



progenitors, multi-dimensional phenomena such as rotation and the development of fluid 



instabilities are required and help to increase the neutrino heating efficiency 



models are also required to aid the understanding of aspherical explosions 37 



33 



- 361] . Such 



III. THE QUARK BAG AND PNJL HYBRID MODELS 



We discuss two different quark matter descriptions for use in astrophysical applications 
as well as their differences in the resultin g E oS. The quark bag model for three flavor quark 



matter is described in detail in the Refs. 



39 



4]J . Bag constants in the range of S 1 / 4 =155- 



165 MeV and corrections from the strong interaction were adopted and a fixed strange quark 
mass of 100 MeV is applied in these simulations which we want to contrast here with first 
results for a PNJL model. A serious drawback of the bag model parametrization employed 
here is that cannot describe the mass of the recently observed 1.97 M neutron star 42]. 
This could be cured by extending the bag model and accounting for leading order QCD 
corrections and diquark condensation 43|, |44j but results in an early onset of quark matter 
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only for considerably stiff nuclear EoS. 

The three flavor quark matter PNJL model is based on the description in Ref. 
with the Polyakov loop extension according to 14j. An additional isoscalar vector meson 
interaction leads to a stiffening of t he q uark matter equation of state and allows to describe 
hybrid stars with a mass of 2 M 45]. The diquark and vector meson coupling is set to 
rip = 1.02 and rjy = 0.25, respectively. Strangeness on the quark side in the PNJL model 
occurs at higher densities than the onset of deconfinement, due to the dynamical quark 
masses involved. The resulting phase diagram is discussed in Ref. p]. 



The Figs. 1(a) and 1(b) compare the phase diagrams for the quark bag model and the 
PNJL model for different proton-to-baryon ratios Y p relevant for supernova matter. There 
are different critical densities for the onset of deconfinement comparing the quark bag and 
PNJL models for equal Y p . E.g. for Y p = 0.3 and T = 0, the critical densities are 0.159 fm _1 
for the bag model in comparison to 0.214 fm _1 for the PNJL model (see the solid lines 



in the Figs. 1(a) and 1(b)). For the bag model, the early onset of deconfinement close 
to normal nuclear matter density illustrates the strong isospin dependency of the critical 
density. For instance, at higher Y p = 0.5 the critical density shifts to 0.321 fm _1 for T = 0. 
In general, small transition densities result from the application of the Gibbs conditions for 
the phase transition with the corresponding extended mixed phase and from the existence of 
the s-quark flavor which results in a larger number of degrees of freedom when compared to 
nuclear matter. For the PNJL model, the transition to quark matter, shown in Fig. 1(b) , is 



based on the Maxwell construction. The critical densities are generally larger and the mixed 
phase is very narrow in density compared to the quark bag model. The PNJL model also 
shows the typical isospin dependency. The critical density reduces at decreasing Y p , which 
is however much weaker than for the quark bag model. Both models differ also significantly 
in the maximum mass of neutron stars, 1.50 M for the bag model with B 1 ^ = 165 MeV 
and 1.97 M Q for the PNJL model. 

Furthermore, the temperature dependence of the critical density is different in both quark 
matter descriptions. This is of particular relevance for astrophysical applications that explore 
a possible quark-hadron phase transition. In the quark bag model, the critical density 
reduces continuously at higher temperatures. It reaches about 0.1 fm _1 at T = 50 MeV for 
Y p = 0.3. For the PNJL model, the critical density rises with increasing temperature for 
any Y p . This behavior is a result of the neglect of any modification of the Polyakov loop 
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potential on the quark density and shou 



such a procedure can be defined, see 



d be improved as soon as reasonable constraints for 



46( for a first step. 



The resulting EoS are illustrated in Figs. 2(a) and 2(b) , showing the pressure-density 



curves for different entropies per baryon and fixed Y p = 0.3. The largely extended mixed 
phase region in the phase diagram for the bag model corresponds to a soft EoS. This is 



illustrated via the pressure-density curves in Fig. 2(a), where the adiabatic index differs 
initially only little between the onset of deconfinement and the pure hadronic phase. To- 
wards the end of the mixed phase where the pure quark phase sets in, the adiabatic index 
is largely reduced for any entropy per baryon and hence the EoS is significantly softer com- 
pared to the hadronic case. In the pure quark phase, the adiabatic index increases again 
significantly, which makes the pure quark phase stiffer than the mixed phase. The sharp 
transition between the mixed and the pure quark phases, gives rise to a strong effect for the 
hydrodynamics evolution along isentropes. 



The PNJL hybrid EoS shows a different pressure behavior (see Fig. 2(b) ) during the 
quark-hadron phase transition. The largest change of the adiabatic index is found at the 
onset of deconfinement. Close to the onset of pure quark matter, the adiabatic index changes 
only little. These two aspects are a general feature of the chosen construction of the PNJL 
hybrid EoS, for which a smooth transition can be expected in dynamical simulations. 



IV. CORE COLLAPSE SUPERNOVA EVOLUTION IN THE PHASE DIAGRAM 



The central temperature and density evolution of the 15 M core collapse supernova 



simulation is shown via the dashed lines in the Figs. 1(a) and 1(b), for the first second 
post bounce, during which the explosion is expected to take place for such massive Fe-core 
progenitors [38]. This phase of the supernova is determined by mass accretion (~ 0.1 M /s) 
onto the central PNS, which consequently contracts on a timescale between 100 ms up to sec- 
onds. During the PNS contraction, the central density and temperature rise continuously 2 . 
During this evolution, the central electron fraction reduced from ~ 0.3 to Y p ~ 0.25. 



The simulation is based on a pure hadronic description of matter 



2j. The current super- 



nova models explore only hadronic EoS. It is shown here to illustrate the possibility and the 



2 The highest temperatures ~ 30 — 60 MeV are not obtained at the center of the PNS but at slightly lower 
densities; it corresponds to the region where the bounce shock formed initially. 
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conditions relevant in order to obtain quark matter at supernova cores. 

The differences between the two hybrid EoS become clear. Using the quark bag model, 
it is possible to obtain quark matter. The early post bounce evolution of the central mass 
trajectories considered for the 15 M progenitor model, enter deeply into the mixed phase. 
The same holds true for the non-central part of the PNS, where up to 0.5 M Q can reach the 
mixed phase within the first 500 ms post bounce evolution (depending on the progenitor 



model 



4=11]) . Within the PNJL model, the central densities obtained for this particular 15 M 



core collapse supernova simulation are not sufficiently high enough to enter the mixed phase. 
It would require an additional rise of the central density by a factor of two or more. For 
such quark-hadron hybrid EoS quark matter will not be reached in core collapse supernovae 
of massive stars up to 15 M Q as discussed here, during the expected explosion phase of 
0.5-1 seconds post bounce. 

In order to investigate the EoS softening effects from the presence of quark matter in the 
PNS interior and possible consequences for the dynamical evolution (and possible observa- 
tions), we apply the quark bag model hybrid EoS (using B 1 ^ = 155 MeV, as = 0.3) to a core 
collapse supernova evolution (again the 15 M model). The largely reduced adiabatic index 
at the end of the mixed phase causes a gravitational collapse of the PNS, spatially separated 
into central sub-sonic collapse and outer super-sonic collapse. During the collapse density 
and temperature rise (see Fig. [3]), which in turn favors pure quark matter over hadronic mat- 
ter. In the pure quark phase, where the adiabatic index is increased, the collapse halts and 
a strong hydrodynamic shock front forms. This was first observed in the context of a first 



order deconfinement phase transition in ref. y. 
multi-peaked neutrino signal form SN1987A 



, in relation to the statistically insignificant 



refs. 



50 



491 ]. and later discussed in more detail in the 



51] . The shock wave appears initially as a pure accretion front, which is illustrated 



in Fig. [3] at 310.4668 ms post bounce. The shock propagates outwards in radius towards 
the PNS surface, driven by the thermal pressure of the deconfined quarks and neutrino 
heating behind the shock. The latter aspect is related to local heating rates due to electron 
(anti)neutrino absorptions, which increase at the quark-hadron phase boundary where the 
infalling hadronic material converts into quark matter and density and temperature rise by 
several orders of magnitude. At the PNS surface, where the density decreases over several 
orders of magnitude, the accretion front accelerates and positive velocities are obtained (see 
Fig. [3] at 310.5135 ms post bounce). This moment determines the onset of explosion, even 
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in core collapse supernova models where otherwise (i.e. without the quark-hadron phase 
transition) no explosions could have been obtained. The shock expands and continues to 
accelerate where matter velocities on the order of 10 5 km/s are obtained. 

The shock propagation across the neutrinospheres releases an additional burst of neutri- 
nos, illustrated in Fig. HI This second burst rises in all neutrino flavors. It differs from the 
reference case discussed in ref. [411 ]. where the second burst was dominated by u e and (z^/ T , 
v^/r) due to the lower density of the PNS envelope of the less massive progenitor chosen 
and the different quark bag EoS parameters. The rise of the burst is due to the presence of 
numerous electron-positron pairs that allow for the production of v e and v e via the charged 
current reactions as well as (zv/ r , zV/ r ) v i a pair processes. It differs from the deleptonization 
burst at core bounce, which is only in v e emitted via a large number of electron captures. 
This second burst, if occurring, is a strong observational indication for the reconfiguration 
of the high density domain in core collapse supernova evolution. 

Coming back to the PNJL model in core collapse supernova simulations, preliminary 
results of long-term non-exploding models in the progenitor mass range between 15-25 M Q , 
show a smooth transition to quark matter at the PNS interior where only the mixed phase 
could be reached. Very little mass (~ 0.05 M ) is converted into the mixed phase and the 
quark fraction rises only on timescales on the order of seconds, which corresponds to the 
PNS contraction timescale induced via mass accretion. However, massive star explosions 
are expected to take place during the first 0.5-1 seconds post bounce. Hence, a direct 
correlation to the explosion mechanism of massive stars as well as direct observables that 
relate to the explosion phase, cannot be expected for such a hybrid EoS. However, this hybrid 
EoS model has another interesting feature. Since the PNJL model possesses a critical point, 
one may study whether in the collapse of very massive progenitors with subsequent black 
hole formation the trajectory of the collapse in the QCD phase diagram sweeps the critical 
point before horizon crossing [52]. In that case a characteristic modification of possible 
observables from the deconfinement transition is expected, at variance to trajectories which 
do not pass the critical point. The realistic modeling of such processes requires further 
development of a phase transition construction which does not remove the critical endpoint. 



10 



V. SUMMARY 

We discussed the possibility to reach the critical conditions for the onset of deconfinement 
in core collapse supernovae, comparing two different quark-hadron hybrid EoS. Based on 
the bag model for strange quark matter and choices of low bag constants, without violating 
constraints from e.g. heavy-ion collision experiments, the low critical densities as well as the 
strong isospin and temperature dependencies allow for the quark-hadron phase transition to 
take place in core collapse supernova simulations within the first 500 ms post bounce. For a 
particular choice of parameters, we illustrated the dynamical evolution. It was determined 
by an adiabatic collapse due to the soft EoS in the extended mixed phase modeled by 
applying Gibbs conditions, and the formation of a strong hydrodynamic shock wave. It 
triggered the explosion even in models where otherwise explosions could not be obtained. 
Furthermore, the millisecond neutrino burst released from the shock propagation across the 
neutrinospheres may become observable for a future Galactic event if quark matter occurs 



53|. The future observation of such a multi-peaked neutrino signal may reveal information 
about the nuclear EoS at high densities and temperatures that cannot be reached in heavy- 
ion collision experiments at present. For the PNJL model parametrization presented here, 
the higher critical density and the narrow mixed phase due to the Maxwell construction, 
together with the rising critical density for increasing temperatures, prevent the quark- 
hadron phase transition to occur in core collapse supernova simulations during the post 
bounce accretion phase. For such an EoS, a direct correlation between the phase transition 
to quark matter and the explosion mechanism cannot be expected. However, the conditions 
for the appearance of quark matter for the PNJL model may allow to obtain a signal from 
core collapse supernovae during the PNS deleptonization phase, on timescales on the order 
of several seconds after the onset of, e.g., initially neutrino-driven explosions. 

The two classes of hybrid EoS discussed here allow for very different types of investigations 
in simulations of core collapse supernovae. Hybrid EoS based on the class of bag models for 
quark matter give access to bulk properties of quark-hadron matter and allow to model EoS 
which undergo a quark-hadron phase transition during the supernova collapse accompanied 
by a second neutrino burst. However, the bag model cannot capture the aspect of chiral 
symmetry restoration and possible effects of the existence of a critical point. Moreover, 
it remains to be shown whether the 2 M mass constraint can be explained in agreement 
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with the deconfinement scenario for this class of models. The hybrid models based on a 
PNJL-type approach to quark matter would allow to study the question of the existence 
of critical points in the QCD phase diagram, e.g., during black hole formation where high 
temperatures on the order of several 100 MeV are obtained. A systematic study of possible 
observables from black hole formation using PNJL models that include a critical endpoint, 
may provide the astrophysical analogue of the energy scan programs in heavy-ion collisions. 
However, the 2 M Q mass constraint is obtained at the price of a too high critical density 
for a deconfinement transition during collapse. This statement is not the final word. It is 
possible, e.g., via a density dependence of the gluon thermodynamics which is absent in the 
present model, to allow for a simultaneous description of deconfinement in the supernova 
interior and high mass (proto) neutron stars. 
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(a) Quark bag model (B 1 / 4 = 165 MeV), for 
Y p — 0.1 (solid lines), Y„ — 0.3 (dash-dotted lines) 
and Y p = 0.5 (dotted lines) 
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(b) PNJL model (i}n = 1-02, r/ v = 0.25), for 
Y p = 0.2 (solid lines), Y„ — 0.3 (dash-dotted lines) 
and Y p = 0.5 (dotted lines) 



FIG. 1: QCD phase diagrams for the two hybrid EoS under discussion (onset of quark matter: 
thin lines, pure quark phase: thick lines) for selected proton-to-baryon ratios Y p that are indicated 
via the text arrows in the graphs. The dashed line shows the evolution of the central temperature 



and density for a representative core collapse supernova of a 15 M Q progenitor model [47], based 
on a hadronic EoS [2] for the first second post bounce. Note the kink between 0.12-0.15 fm -1 , it is 
related to the core bounce after which the central density and temperature decrease slightly before 
they continue to increase during the later compression. The wide stretched mixed phase region for 
the bag model in graph (a) results from the applied Gibbs construction while the in comparison 
narrow region for the PNJL model in graph (b) is a result of the Maxwell construction. 
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(a) Quark bag model (B 1 / 4 = 165 MeV). 
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(b) PNJL model {r) D = 1.02, r] V = 0.25). 



FIG. 2: Total pressure density curves for different entropies per baryon and fixed Y p = 0.3, based 
on the two hybrid EoS. 
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FIG. 3: Radial profiles of velocity, density and temperature for a 15 M Q progenitor model at 
selected post bounce times during the QCD phase transition induced PNS collapse, using the bag 
model with B 1 / 4 = 155 MeV and a s = 0.3. 
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FIGURE CAPTIONS 

1. QCD phase diagrams for the two hybrid EoS under discussion (onset of quark matter: 
thin lines, pure quark phase: thick lines) for selected proton-to-baryon ratios Y p that 
are indicated via the text arrows in the graphs. The dashed line shows the evolution 
of the central temperature and density for a representative core collapse supernova 
of a 15 M progenitor model 47J], based on a hadronic EoS |2j| for the first second 
post bounce. Note the kink between 0.12-0.15 fm -1 , it is related to the core bounce 
after which the central density and temperature decrease slightly before they continue 
to increase during the later compression. The wide stretched mixed phase region for 
the bag model in graph (a) results from the applied Gibbs construction while the in 
comparison narrow region for the PNJL model in graph (b) is a result of the Maxwell 
construction. 

2. Total pressure density curves for different entropies per baryon and fixed Y p = 0.3, 
based on the two hybrid EoS. 

3. Radial profiles of velocity, density and temperature for a 15 M Q progenitor model at 
selected post bounce times during the QCD phase transition induced PNS collapse, 
using the bag model with B 1 / 4 = 155 MeV and as = 0.3. 

4. Evolution of the neutrino luminosities and mean energies for a 15 M model, including 
the quark hadron phase transition using the quark bag model (S 1//4 = 155 MeV, as = 



0.3) |40l . |41| . Due to the phase transition, the PNS becomes gravitational unstable 
and collapses. The critical conditions for the gravitational collapse (T = 17.67 MeV, 
p = 5.511 x 10 14 g/cm 3 , Y p = 0.197) are obtained at about 308 ms after bounce for this 
particular model. The second burst is released when the second shock, that forms due 
to the quark hadron phase transition in the PNS interior, crosses the neutrinospheres. 
The second burst occurs in all neutrino flavors, while the deleptonization burst at core 
bounce appears only in v e . This signature of deconfinement is absent in models that 



explore standard hadronic EoSs (e.g. from ref. [2[). 
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FIG. 4: Evolution of the neutrino luminosities and mean energies for a 15 M model, including the 



40, 



quark hadron phase transition using the quark bag model (Z? 1//4 = 155 MeV, as = 0.3) 
Due to the phase transition, the PNS becomes gravitational unstable and collapses. The critical 
conditions for the gravitational collapse (T = 17.67 MeV, p = 5.511 x 10 14 g/cm 3 , Y p = 0.197) are 
obtained at about 308 ms after bounce for this particular model. The second burst is released when 
the second shock, that forms due to the quark hadron phase transition in the PNS interior, crosses 
the neutrinospheres. The second burst occurs in all neutrino flavors, while the deleptonization 
burst at core bounce appears only in v e . This signature of deconfinement is absent in models that 
explore standard hadronic EoSs (e.g. from ref. 2j). 



